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ABSTRACT 

The 105-ms X-ray pulsar J1852+0040 is the central compact object (CCO) in supernova remnant 
Kes 79. We report a sensitive upper limit on its radio flux density of 12 /zJy at 2 GHz using the NRAO 
Green Bank Telescope. Timing using the Newton X-Ray Multi-Mirror Mission (XMM-Newton) and 
the Chandra X-ray Observatory over a 2.4 yr span reveals no significant change in its spin period. 
The 2cr upper limit on the period derivative leads, in the dipole spin-down formalism, to an energy 
loss rate E < 7 x 10 33 ergs s _1 , surface magnetic field strength B p < 1.5 x 10 11 G, and characteristic 
age r c = P/2P > 8 Myr. This value of r c exceeds the age of the SNR by 3 orders of magnitude, 
implying that the pulsar was born spinning at its current period. However, the X-ray luminosity of 
PSR J1852+0040, L hoi w 3 x 10 33 (d/7.1 kpc) 2 ergs s" 1 , is a large fraction of E, which challenges the 
rotation-powered assumption. Instead, its high blackbody temperature /cTbb = 0.46 ±0.04 keV, small 
blackbody radius Rbb ~ 0.8 km, and large pulsed fraction f p w 80%, may be evidence of accretion 
onto a polar cap, possibly from a fallback disk made of supernova debris. If B p < 10 10 G, an accretion 
disk can penetrate the light cylinder and interact with the magnetosphere while resulting torques on 
the neutron star remain within the observed limits. A weak B-fteld is also inferred in another CCO, 
the 424-ms pulsar IE 1207.4—5209, from its steady spin and soft X-ray absorption lines. We propose 
this origin of radio-quiet CCOs: the magnetic field, derived from a turbulent dynamo, is weaker if 
the NS is formed spinning slowly, which enables it to accrete SN debris. Accretion excludes neutron 
stars born with both B p < 10 11 G and P > 0.1 s from radio pulsar surveys, where B p < 10 11 G is not 
encountered except among very old (r c > 40 Myr) or recycled pulsars. Finally, such a CCO, if born 
in SN 1987A, could explain the non-detection of a pulsar there. 

Subject headings: ISM: individual (Kesteven 79, SN 1987A) — pulsars: individual (IE 1207.4-5209, 
CXOU J185238.6+004020, PSR J1852+0040) — stars: neutron 



1. INTRODUCTION 

A compact X-ray source, CXOU J185238.6+004020, 
was found in the center of the supernova remnant (SNR) 
Kes 79 by Seward et al. (2003). Using XMM-Newton 
in 2004 October, we discovered 105-ms X-ray pulsations 
from this CCO (Gotthelf et al. 2005, Paper 1), also 
named PSR J1852+0040. Its pulsed fraction was as 
high as 86% and it had an apparently thermal X-ray 
spectrum. From the discovery observations, only an up- 
per limit on its period derivative could be determined, 
P < 7 x 10~ 14 s s _1 , which left the energetics of the 
NS and the mechanism of its X-ray emission ambigu- 
ous. The two most plausible models explored in Paper 
1, a rotation-powered pulsar, and accretion from fallback 
material, each encountered some difficulties. 

Here, we present three new observations, two from 
XMM-Newton and one from Chandra, that over a 2.4 
year span refine the spectral and timing properties of 
PSR J1852+0040. In addition, we report the nega- 
tive result of a search for radio pulsations in a deep 
GBT pointing. From the assembled X-ray timing of 
PSR J1852+0040, a sensitive new upper limit is obtained 
on its period derivative, suggesting that the underlying 
cause of its unusual properties is a weak magnetic field, 
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B p < 1.5 x 10 11 G. Unlike canonical radio pulsars with 
Bp ~ 10 12 ~ 13 G, a young NS that is weakly magnetized 
and spinning slowly at birth is able to accrete from ma- 
terial that is thought to remain following a supernova 
explosion. We discuss this as a model for the proper- 
ties of radio-quiet CCOs, including the similar pulsar 
IE 1207.4-5209. 

2. XMM-Newton AND Chandra OBSERVATIONS 

The two observations previously reported in Paper 1 
were obtained with XMM-Newton on 2004 October 18 
and 23. Two more observations, reported here, were 
made on 2006 October 8 and 2007 March 20 using the 
same instrument modes, and exposure times of 30 ks 
each. The pn CCD of the European Photon Imag- 
ing Camera (EPIC; Turner et al. 2003), was operated 
in "small window" mode with a 4f3 x 4'3 field-of-view 
(FOV) and 5.7 ms time resolution to resolve the pulsa- 
tions of PSR J1852+0040. The EPIC MOS CCDs were 
exposed in "full frame" mode with a 30' diameter FOV 
and time resolution of 2.7 s. For all three cameras the 
medium density filter was used. We processed all four 
observations with Science Analysis System version xmm- 
sas_20060628_1801-7.0.0, using the same techniques as in 
Paper 1. 

A 32 ks Chandra observation was performed on 2006 
November 23 using the Advanced Camera for Imag- 
ing and Spectroscopy (ACIS) operated in continuous- 
clocking (CC) mode to provide a time resolution of 
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Fig. 1. — Geometry of the 2006 November 23 Chandra CC-mode 
observation of PSR J1852+0040, showing the outline of the S3 
CCD and the columns used for source and background extraction, 
overlaid on the 2001 July 31 ACIS-I observation of Seward ct al. 
(2003). (This broadband image of Kes 79 is square-root scaled to 
emphasize faint emission features and background.) The counts 
were summed along the length of the columns, oriented at posi- 
tion angle 241°.5 east of north. Both spectral and timing data for 
the pulsar were extracted from the six central columns. The back- 
ground spectrum was obtained from nine neighboring columns on 
either side (flanking lines), separated from the source by a two- 
column gap. 



2.85 ms. The back-illuminated ACIS-S3 CCD was used, 
and the Scientific Instrument Module (SIM) was offset in 
the — Z direction to place the pulsar 0'.85 from the edge 
of the CCD to reduce contamination from the thermal 
remnant. To achieve the fast timing in CC mode, one 
spatial dimension of the CCD image is sacrificed by in- 
tegrating along the column direction during each CCD 
readout cycle. The geometry of the observation is illus- 
trated in Figure 1, overlaid on an image of Kes 79 previ- 
ously obtained with ACIS-I in timed exposure mode on 
2001 July 31 (Seward et al. 2003). We also reanalyzed 
the 2001 observation for the spectral properties of the 
pulsar. During both observations the CCD background 
count rate was stable and no time filtering was necessary. 
All photon arrival times were adjusted by the standard 
processing to account for the spacecraft dithering, and 
SIM offset in the case of the CC-mode observation. Re- 
duction and analysis of the Chandra data are all based on 
the standard software package CIAO (v3.4) and CALDB 
(v3.3). 

2.1. Timing Analysis 

Photon arrival times were converted to the solar system 
barycenter using the Chandra derived source coordinates 
from Paper 1, R.A. = 18 h 52 m 38!57, decl. = +00°40'19"8 
(J2000.0), and signal strength was maximized by select- 
ing the 1 — 5 keV energy band. From the XMM-Newton 
observation of 2006 October 8, a Z\ test (Bucchcri et al. 
1983) localizes the period to P = 104.912593(20) ms with 
a peak statistic of Z\ = 88. The resulting light curve, 
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Fig. 2.— Folded light curves of PSR J1852+0040 in the 1 - 5 
keV band using XMM-Newton EPIC pn on 2006 October 8 (top), 
Chandra ACIS-S3 in CC mode on 2006 November 23 (middle), 
and XMM-Newton EPIC pn on 2007 March 20 (bottom). The 
estimated background level is indicated by the dashed line. Phase 
zero is arbitrary. 



which has pulsed fraction f p = 77% ± 15% after account- 
ing for SNR background, is shown in Figure 2. Here, we 
define the pulsed fraction as f p = A(pulsed)/A(total), 
where we choose the minimum of the folded light curve as 
the unpulsed level. The observation on 2007 March 20 
yielded consistent parameters, P — 104.912609(19) ms 
with Z\ = 83 and f p = 71% ± 15%. These new values of 
f p are intermediate between the 2004 XMM-Newton ob- 
servations, which yielded f p = 86%±16% and 61%±16%, 
respectively, weakening the marginal evidence of variabil- 
ity of f p that was noted in Paper 1. Similar analysis of 
the Chandra observation of 2006 Novemeber 23 yields a 
formally identical period, P ~ 104.912612(19) ms with 
Z\ = 94, and f p consistent with 100%. The pulsed frac- 
tion is less reliable in CC-mode data due to systematics 
of estimating the contamination from SNR background. 
Table 1 is a summary of all of the available timing data 
on PSR J1852+0040'. 

While the new observations extend the sampled time 
span by a factor of 150 since Paper 1, there is still no 
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TABLE 1 
X-ray Timing of PSR J1852+0040 



Mission 


Date 


Epoch 


Duration 


Ratc a 


Period b 


f P c 




(UT) 


(MJD) 


(s) 


(s-i) 


(ms) 


(%) 


XMM 


2004 Oct 18 


53,296.001 


30,587 


0.041(2) 


104.912643(18) 


86(16) 


XMM 


2004 Oct 23 


53,301.985 


30,515 


0.046(2) 


104.912600(27) 


61(16) 


XMM 


2006 Oct 08 


54,016.245 


30,243 


0.045(2) 


104.912593(20) 


77(15) 


Chandra 


2006 Nov 23 


54,062.257 


32,165 


0.022(2) 


104.912612(19) 


100(16) 


XMM 


2007 Mar 20 


54,179.878 


30,506 


0.043(2) 


104.912609(19) 


71(15) 



a For XMM-Newton, background and dead-time corrected count rate in a 15" radius aperture. Statistical (V~N) uncertainty in the last 
digit is given in parentheses. 

b Period derived from a Z\ test. Period uncertainty is la computed by the Monte Carlo method described in Gotthelf et al. (1999). 
c Pulsed fraction defined as f p = A r (pulscd)/7V (total). Unpulsed level set at lowest point of 20-bin folded lightcurve. 



TABLE 2 

X-ray Spectral Fits for PSR J1852+0040 



Parameter 


2001 Jul 31 

Chandra 


2004 Oct 18 
XMM 


2004 Oct 23 
XMM 


2006 Oct 08 
XMM 


2006 Nov23 
Chandra 


2007 Mar 20 
XMM 


Blackbody Model 


N-a (10** cm"*) 
kT BB (keV) 
Rbb (km) 
Abb (10 10 cm 2 ) 
Fbb(10- 13 cgs) a 
i BB ,boi(10 33 cgs) b 
X» 


1.1 ±0.3 
0.50 ± 0.04 
0.62 ± 0.05 

4.8 ±0.7 
2.0 ±0.3 

3.2 ± 0.7 
0.71(18) 


1.4 ±0.3 
0.44 ± 0.04 
0.92 ± 0.09 
10.7 ± 2.1 
1.9 ±0.3 
4.3 ± 1.1 
0.85(68) 


1.4 ± 0.3 
0.45 ± 0.03 
0.79 ± 0.09 
7.8 ± 1.7 
2.0 ± 0.3 
3.3 ± 1.0 
0.67(70) 


1.4 ± 0.3 
0.46 ± 0.04 
0.74 ±0.08 
6.8 ± 1.5 
2.0 ± 0.3 
3.2 ± 1.0 
0.86(69) 


I. 6 ±0.3 
0.41 ± 0.04 
0.96 ±0.04 

II. 5 ±4.4 
1.7±0.4 
3.5 ± 1.7 
0.88(22) 


1.3 ±0.3 
0.46 ± 0.04 
0.69 ± 0.09 
6.0 ± 1.5 
1.9 ±0.3 
2.9 ±1.1 
0.67(69) 



Power-law Model 



JV H (10 22 cm" 2 ) 


2 6 +a7 


o 9 +0.6 
•^-0.6 


O 1+0.6 

• s - 1 -0.5 


q 9+O.6 


o 1+0.8 
J - 4 -0.6 


q 9+0.7 


r 


4 1+0-7 


4 9+0-7 


4 O+06 


4 8 +0 ' 6 


c q+0.8 


4 n+0.6 


Fp L (10- 13 cgs) a 
L PL (10 33 cgs) c 


2.0 ±0.6 


1.9 ±0.43 


2.0 ±0.4 


2.1 ±0.5 


1.8 ±0.4 


1.9 ±0.4 


1.2 ±0.4 


1.1 ±0.2 


1.2 ±0.2 


1.2 ±0.3 


1.1 ± 0.3 


1.2 ±0.3 




0.94(18) 


0.82(68) 


0.68(70) 


0.78(69) 


0.81(22) 


0.71(69) 



Note. — Uncertainties are 68% confidence intervals for two interesting parameters. Luminosities are computed for d ■ 

a Absorbed flux in the 1 — 5 keV band in units of ergs cm -2 s -1 . 
b Unabsorbed, bolometric blackbody luminosity in units of ergs s -1 . 

c Absorbcd luminosity in the 1 — 5 keV band in units of ergs s _1 . This corrects an error in Table 2 of Paper 1. 



7.1 kpc. 



significant detection of a period derivative. A x 2 fit to 
the five measurements of P yields the formally negative 
< P >= (-3.4 ± 2.7) x f(T 16 s s" 1 that is, however, 
consistent with zero at the w la level. We are unable to 
link any of the observations with a definite cycle count 
in order to improve on this measurement. If we adopt 
the 2a upper limit P < 2.0 x 10~ 16 , the dipole spin- 
down formalism for an isolated pulsar implies an energy 
loss rate E = -Ifltl = 4ir 2 IP/P 3 < 7 x 10 33 ergs s"\ 

surface magnetic field strength B p = 3.2 x 10 19 \/ PP < 
1 .5 x fO 11 G, and characteristic age r c = P/2P > 8 Myr. 
These derived quantities pose problems for a rotation- 
powered X-ray source, as discussed in §4.1. 

2.2. Spectral Analysis 

Following the analysis method in Paper 1 for each 
XMM-Newton observation, spectra from the two EPIC 
MOS cameras (MOS1, MOS2) were co-added. The 
EPIC MOS and EPIC pn spectra were then fitted simul- 



taneously using the XSPEC (vl2.21) package. The results 
of spectal fits using either a blackbody or an absorbed 
power-law model are presented in Table 2. An acceptable 
X 2 statistic is obtained using either model. However, the 
blackbody is preferred over the power law based on its 
derived column density of N B = (1.4 ± 0.3) x 10 22 cm -2 , 
which is consistent with that found for the SNR (iVu ~ 
1.6 x 10 22 cm" 2 , Sun ct al. 2004). The fitted column den- 
sity for the power-law model, A^h = 3.2lgg x 10 22 cm -2 , 
is significantly larger than the integrated 21 cm Galactic 
value of 2 x 10 22 cm -2 in this direction (Dickey & Lock- 
man 1990). The best-fit blackbody model yields tem- 
perature kT BB = 0.46 ± 0.04 keV (see Fig. 3). Adopt- 
ing a distance d = 7.1 kpc from Case & Bhattacharya 
(1998), the bolometric blackbody luminosity iBB.boi ~ 
3.0 x 10 33 ergs s _1 , corresponding to blackbody area 
A BB w 7xl0 10 (rf/7.1 kpc) 2 cm 2 or w 0.4% of the NS sur- 
face. These results are consistent with the spectral pa- 
rameters from the previous XMM-Newton observations 
in 2004 October, described in Paper 1 and summarized 
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Fig. 3.— X-ray spectra of PSR J1852+0040 in Kes 79. Top: 
The XMM-Newton spectrum on 2006 October 8. Data from the 
EPIC pn {black, upper) and EPIC MOS (red, lower) with the best- 
fit blackbody model using parameters in Table 2. The residuals 
from the fits are in units of a. Middle: Spectrum on 2006 November 
23 from the Chandra ACIS-S3 in CC mode. Bottom: The XMM- 
Newton spectrum on 2007 March 20. [See the electronic version of 
the Journal for a color version of this figure.] 



here in Table 2, indicating that the flux and spectrum 
have remained steady. 

The Chandra spectra from the 2001 and 2006 observa- 
tions were each prepared by grouping a minimum of 40 
counts per channel, and fitted using the XSPEC software. 



The standard spectral response matrices were generated 
for the target location on the CCD and all spectra were 
corrected for the effects of charge-transfer inefficiency; 
however, the spectral gain is not calibrated well in CC- 
mode. The resulting spectral parameters from the Chan- 
dra observations, listed in Table 2, are similar to the 
ones from XMM-Newton. The blackbody fits to the new 
XMM-Newton and Chandra spectra are shown in Fig- 
ure 3. 

3. SEARCH FOR A RADIO PULSAR 

On 2002 October 29 we observed the location of 
CXOU J185238.6+004020 with the ATNF Parkes tele- 
scope in NSW, Australia. At that time the period of 
the pulsar was not yet known, so we employed standard 
pulsar searching techniques. We used the central beam 
of the multibeam receiver at a frequency of 1374 MHz, 
with 96 channels across a bandwidth of 288 MHz in each 
of two polarizations. The integration time was 6.8 hr, 
during which total-power samples were obtained every 
0.73 ms and recorded for off-line analysis. We searched 
the dispersion measure range 0-2500 cm~ 3 pc (twice the 
maximum Galactic DM predicted for this line of sight by 
the Cordes & Lazio 2000 electron density model; for the 
distance of 7.1 kpc, the predicted DM is 440cm _3 pc). 
The search followed closely that described in more detail 
in Camilo et al. (2006). No promising candidate periods 
were identified. 

Following the X-ray discovery of the pulsar period, we 
made a deeper search using the GBT. We did this on 
2005 June 10, observing the pulsar for 11.1 hr with the 
Spigot spectrometer (Kaplan et al. 2005), sampling 768 
frequency channels across a bandwidth of 600 MHz cen- 
tered on 1950 MHz every 81.92 fts. In addition to per- 
forming a standard blind search, we also folded the data 
at a number of DMs for periods close to the X-ray value 
to increase the sensitivity, but again found no candidates. 
We used the PRESTO software package to analyze these 
data (Ransom 2001; Ransom et al. 2002). 

Applying the standard modification to the radiometer 
equation, assuming a quasi-sinusoidal pulse shape (as in 
X-rays; see Fig. 2), and accounting for a sky tempera- 
ture at this location of 5 K, we were sensitive to a pulsar 
with period 105 ms having a period-averaged flux den- 
sity at 2 GHz > 12 /iJy. Converting this to the more 
common 1.4 GHz pulsar search frequency, using a typi- 
cal spectral index of -1.6 (Lorimer et al. 1995), results 
in 5i.4 < 0.02 mJy. For a distance of ~ 7 kpc, this corre- 
sponds to a pseudo-luminosity limit of L1.4 = S1.40P < 
1 mJy kpc 2 . Only one young radio pulsar has a smaller 
radio luminosity than this (the pulsar in 3C58; Camilo 
et al. 2002). 

Based on these results, it is therefore somewhat un- 
likely that PSR J1852+0040 is presently a radio pulsar 
beaming toward the Earth. More likely, it is either beam- 
ing away (observed radio pulsars with similar periods, or 
ages < lOkyr, have beaming fractions of ~ 0.3; see Tau- 
ris & Manchester 1998), or not an active radio pulsar at 
all. 
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4. DISCUSSION 
4.1. Rotation-Powered Pulsar? 

Potential obstacles to the interpretation of PSR 
J 1852+0040 as a rotation-powered pulsar were identified 
in Paper 1. The newly reduced upper limit on its P, by 
more than 2 orders of magnitude, turns those difficul- 
ties into major objections. First, the X-ray luminosity of 
PSR J1852+0040, L bo i « 3 x 10 33 (d/7.1 kpc) 2 ergs s" 1 , 
is comparable to the 2a upper limit on its spin-down 
power, E < 7 x 10 33 ergs s _1 . Thus, it may be difficult 
to power the X-ray source with spin-down energy. Since 
the X-ray spectrum is consistent with thermal emission, 
what about residual cooling of the neutron star? While 
the X-ray luminosity of PSR J1852+0040 is consistent 
with minimal NS cooling curves for an age of 10 3-4 yr 
(Page et al. 2004), its blackbody temperature implies 
an emitting area that is just « 0.4% of the NS surface. 
Even taking into account that, in a real NS atmosphere, 
a blackbody fit overestimates the effective temperature 
and underestimates the area, the highly modulated X- 
ray pulse must come from a small hot spot. Although 
anisotropic conduction in a magnetized NS can enforce 
temperature gradients on the surface, such models (Gep- 
pert et al. 2004; Perez-Azorin et al. 2006) do not show 
spots as small and as hot as observed here. Instead, an 
external source of localized heating may be needed. 

The open-field-line polar cap is the most likely target 
area for external heating of the NS surface, with the en- 
ergy ultimately drawn from rotation. The canonical area 
of the polar cap is A pc = 2ir 2 R 3 /Pc kIx 10 10 cm 2 ; this 
is ~ 14% of the area implied by the blackbody fit to the 
X-ray spectrum. Among theories of polar-cap heating, 
a maximum thermal luminosity is predicted by Wang et 
al. (1998) from the outer-gap model for 7-ray pulsars. 
In this analysis, the X-ray luminosity of the hot polar 
cap is limited by the Goldrcich- Julian e^ current flow of 
N « 2 x 10 31 (P/0.105 s)~ 2 (Bp/10 11 G) s" 1 depositing 
an average energy per particle Ef w 4.3 ergs. The max- 
imum luminosity is Lboi ~ fEfNo < 4 x 10 31 ergs s _1 . 
Here, the fraction / of the current reaching the surface 
is set to its estimated maximum value of |. Polar-cap 
heating models of Harding & Muslimov (2001, 2002) pre- 
dict even less X-ray luminosity than Wang et al. (1998). 
Therefore, current theories fall short of predicting the 
temperature and luminosity of the X-ray emission from 
PSR J1852+0040 in the context of a spinning-down neu- 
tron star. 

After the unexplained X-ray properties, a lesser prob- 
lem for PSR J1852+0040 as an isolated pulsar is its char- 
acteristic age r c = P/2P > 8 Myr, which exceeds by 3 
orders of magnitude the SNR age, estimated dynamically 
as 5.4 - 7.5 kyr (Sun et al. 2004). This would require the 
pulsar to have been born at its current period to be as- 
sociated with the SNR. While not fundamentally contra- 
dicting the rotation-powered hypothesis, this result is a 
definite example in support of a recent population anal- 
ysis that favors a wide distribution of radio pulsar birth 
periods (Faucher-Giguere & Kaspi 2006), even though 
PSR J1852+0040 is not itself a radio pulsar. Further- 
more, as magnetic field is generated by a turbulent dy- 
namo whose strength depends on the rotation rate of 
the proto-neutron star (Thompson & Duncan 1993), it is 



natural that long-period pulsars would have the weaker 
£>-ficlds at birth, and the model of Bonanno et al. (2006) 
shows this. Bonanno et al. (2006) also find that a slowly 
rotating NS should have its £?-field confined to small- 
scale regions, with the global dipole that is responsible 
for spin-down absent. 

Finally, we note that there are no young radio pul- 
sars known with B p < 10 11 G. All such weakly magne- 
tized radio pulsars are either very old, with r c > 40 Myr 
(Manchester et al. 2005) 3 , or are recycled, the result of 
accretion from a binary companion. Faucher-Giguere & 
Kaspi (2006) conclude that there is no evidence for mag- 
netic field decay during the radio emitting life of a pulsar, 
so the observered distribution of B p should resemble its 
birth distribution. Among the B p distribution of ordi- 
nary (not recycled) radio pulsars, PSR J1852+0040 is 
certainly in the bottom 10%, and maybe lower. PSR 
J 1852+0040 may become a radio pulsar in the far fu- 
ture, but its weak B-field may be preventing radio emis- 
sion now by permitting accretion. While there is no 
evidence of a close binary companion in the timing of 
PSR J1852+0040 (see Paper 1), it is possible that low- 
level accretion of SN debris prevents it and other radio- 
quiet CCOs from becoming radio pulsars for thousands 
or even millions of years. We next consider whether the 
hypothesis of accretion from such a fallback disk (e.g., 
Alpar 2001; Shi & Xu 2003; Ek§i et al. 2005) better ex- 
plains the X-ray timing and spectral properties of PSR 
J1852+0040. 

4.2. Fallback Accretion? 

The fact that P is so small leads us to review the vi- 
ability of accretion of supernova debris through a disk. 
The X-ray luminosity of PSR J1852+0040 can be pow- 
ered by accretion of to « 3 x 10 13 g s _1 , or only « 0.1 
lunar masses of supernova debris over the past 7.5 kyr. 
The main barrier to disk accretion is the speed-of-light 
cylinder, of radius ri = cP/2ir = 5 x 10 8 cm. If an ac- 
cretion disk cannot penetrate the light cylinder, the NS 
cannot interact with the disk, and it behaves as a isolated 
pulsar. But if B p is as small as 10 10 G, accretion at a 

rate M > 10 13 g s _1 can penetrate the light cylinder, 
since the magnetospheric radius is then 

r m = 3.2 x 10 8 n% 7 M^ 11 (JL^j cm, (1) 

where the magnetic moment fi — B p R?/2 10 28 B Pj i 
G cm 3 . If so, the system is in the propeller regime, in 
which matter flung out from r m at a rate M takes angular 
momentum from the NS, causing it to spin down. The 
propeller spin-down rate is estimated as 

P^2Mr 2 m r 1 P (l-^) (2) 

(e.g., Menou et al. 1999; Ek§i et al. 2005). Here / w 
10 45 g cm 2 is the NS moment of inertia, and P cq = 

3.2 i4 f s 7 M" 377 (M/Mq)~ 5 / 7 s is the equilibrium, or min- 
imum period for disk accretion that is reached when r m < 
r co . The corotation radius is r co = [GM(P/2n) 2 } 1 ^ 3 = 
3.7 x 10 7 cm. 

3 ATNF Pulsar Catalogue, version 1.29, 

http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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While we have not measured a significant P in order 
to determine the needed M, we can adopt the 2a upper 
limit, P < 2.0 x 1CP 16 , as an estimate of P. Combining 
equations (1) and (2), we have in the propeller scenario 

-2/7 



2.2 x 10" 



16 ,,8/7^3/7 



^28 



13 



1 



M 

P 



M5 



0.105 s 



(3) 



We must distinguish here between M, the matter ex- 
pelled that is responsible for the torque on the NS, and 
rh w 3 x 10 13 g s -1 , the matter accreted, which is re- 
sponsible for the X-ray emission from the surface, pre- 
sumably at a magnetic pole of the neutron star. For 
the propeller model to be self-consistent, M must be 
> m, which is possible according to equation (3) as long 
as Bp < 10 10 G. Even so, M does not contribute sig- 
nificantly to the X-ray luminosity because of the much 
weaker gravitational potential at the magnetospheric ra- 
dius r rn from which it is expelled. If B p is as small as 
7 x 10 8 G, then PSR J1852+0040 is not a propeller, but a 
"slow rotator," with r m < r co and P oq < P. In this spin- 
up regime, it would be even more difficult to measure P, 
since 

3/7 



P w -1.3xl(T 18 /$ 7 m% 7 



M 



P 



0.105 s 



(4) 



(Ghosh & Lamb 1979). 

These equations also apply during the prior evolution 
of the pulsar and show that, even if the accretion rate 
was orders of magnitude higher in the past, the spin- 
up and spin-down times P/P are much longer than the 
age of the SNR. There has not been enough time for 
propeller accretion to have spun the pulsar down from 
a much smaller P, so even in this model it was born at 
essentially its present P. 

Thus, we regard either propeller spindown or accretion 
in a weak magnetic field as consistent with the X-ray 
spectral and timing properties of PSR J1852+0040. If 
accreting, its emitting region could be larger than the 
open-field-line polar cap because it would be connected 
to magnetic field lines extending to r m , which is < rg. 
Neither equation (3) nor equation (4) requires there to 
be much torque noise at the level that present data can 
measure. While detected flickering is an indicator of ac- 
cretion in X-ray binaries even in quiescence (Rutledge et 
al. 2007), we do not have strong evidence of variability of 
PSR J1852+0040. But it is not clear that the processes 
in binaries can be extrapolated to accretion at such low 
rates from a fossil disk. Finally, because of the large 
distance and interstellar extinction to PSR J1852+0040, 
as well as its crowded optical field, existing optical data 
do not rule out the presence of a fallback accretion disk 
around PSR J1852+0040 (see Paper 1 for details). 

4.3. Comparison with IE 12074- 5209 and Other 

CCOs 

Difficulties in understanding the timing behavior of 
another CCO, the 424-ms pulsar IE 1207.4-5209 (Za- 
vlin et al. 2004; Woods et al. 2006), were resolved re- 
cently by Gotthelf & Halpern (2007), who showed that 



it is a steady rotator with E < 1.5 x 10 ergs s 
and B p < 3.5 x 10 11 G. This makes IE 1207.4-5209 
very similar to PSR J1852+0040 in its physical param- 
eters. Accretion from a fossil disk of supernova debris 
(Alpar 2001; Shi & Xu 2003; Ek§i et al. 2005; Liu et al. 
2006) was one of the theories considered by Zavlin et al. 
(2004) to explain the now defunct timing irregularities 
of IE 1207.4—5209. But accretion may still be needed to 
explain the radio-quiet nature of IE 1207.4—5209, and 
the details of its X-ray spectrum. However, unlike PSR 
J1852+0040, upper limits on optical/IR emission from 
IE 1207.4—5209 are comparable to what is expected from 
a geometrically thin, optically thick disk accreting at the 
rate required to account for its X-ray luminosity (Za- 
vlin et al. 2004; Wang et al. 2007). Therefore, it may 
prove necessary to invoke a radiatively inefficient flow, 
or perhaps even accretion of solid particles (Cordes & 
Shannon 2006), in order to proceed with this theory for 
IE 1207.4-5209. We propose that this may be the first 
phase in the life of those neutron stars born rotating 
slowly, and with weak magnetic fields. 

A unique phenomenon displayed by IE 1207.4—5209 is 
the set of broad absorption lines in its soft X-ray spec- 
trum. They are centered at 0.7 keV and 1.4 keV (San- 
wal et al. 2002; Mereghetti et al. 2002), and possibly at 
2.1 keV and 2.8 keV (Bignami et al. 2003; De Luca et 
al. 2004), although the reality of the two higher-energy 
features has been disputed (Mori et al. 2005). Proposed 
absorption mechanisms include electron cyclotron in a 
weak (8 x 10 10 G) magnetic field (Bignami et al. 2003; 
De Luca et al. 2004), atomic features from singly ion- 
ized helium in a strong (2 x 10 14 G) field (Sanwal et al. 
2002; Pavlov & Bczchastnov 2005), and oxygen or neon 
in a normal (10 12 G) field (Hailey & Mori 2002; Mori & 
Hailey 2006). A detailed critique of these models is be- 
yond the scope of this paper. We only remark that: (1) 
local values of B measured from X-ray spectral features 
may differ from the global dipole B that is measured by 
timing, especially if the dipole component is weak (Bo- 
nanno et al. 2006), and (2) the attractive feature of the 
electron cyclotron model for IE 1207.4—5209, the equal 
energy spacing of its lines, is in harmony with the present 
work on PSR J1852+0040 which, by different arguments, 
requires a weak surface magnetic field. 

Xu et al. (2003) argued that electron cyclotron har- 
monics of roughly equal depth could be produced be- 
cause, even though the oscillator strength of the har- 
monic is much less than that of the fundamental, reso- 
nant cyclotron radiative transfer could equalize the line 
strengths. Liu et al. (2006) also explained how objec- 
tions to the cyclotron model involving the relative oscil- 
lator strengths can be overcome by assuming that the 
fundamental is optically thick while the harmonic is op- 
tically thin, resulting in comparable equivalent widths. 
They also proposed a specific geometry that accounts 
for the strength of the lines, which also accommodates 
the small pulsed fraction observed from IE 1207.4—5209. 
This involves a column rising along the magnetic axis and 
viewed from the side. Finally, Liu et al. (2006) remarked 
that such a column may be more easily maintained in an 
accretion scenario. 

While the lower-quality spectra of PSR J1852+0040 do 
not yet show evidence of absorption features, cyclotron 
lines could be weak if either the column density in which 
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they are formed or the viewing geometry is not favor- 
able, or especially if the magnetic field is so small as 
to shift them below the observable X-ray band. The 
fitted interstellar column density to PSR J1852+0040, 
TVh = 1.4 x 10 22 cm~ 2 , is much larger than that to 
IE 1207.4-5209, N H = (3 - 10) x 10 20 cm" 2 (Mereghetti 
et al. 2002). Therefore, the soft X-ray region in which 
the strongest lines in IE 1207.4—5209 are seen is partly 
suppressed in the case of PSR J1852+0040. Finally, it 
is possible that the X-ray spectrum of PSR J1852+0040 
is not a pure blackbody as fitted, but is affected by un- 
modelled, broad cyclotron absorption at the low-energy 
end. 

The half dozen CCOs (for a review, see Pavlov et al. 
2004) are similar in their X-ray luminosities and tem- 
peratures. Therefore, they may be compatible with 
the model that we propose for PSR J1852+0040 and 
IE 1207.4-5209. The only piece of evidence that sharply 
contradicts such unification is the apparent historical 
outburst of the CCO in the Cas A SNR that was in- 
ferred from infrared light echos detected with the Spitzer 
Space Telescope (Krause et al. 2005). The geometry and 
sharpness of the echos implies that beamed emission from 
the central source flared for less than a few weeks about 
50 years ago, and the radiation reprocessed into thermal 
emission from dust require an equivalent isotropic flare 
of ~ 2 x 10 46 ergs. Because such energetic outbursts 
are seen on similar timescales from soft gamma-ray re- 
peaters, Krause et al. (2005) proposed that the Cas A 
CCO is a quiescent magnetar, at the other extreme of 
magnetic field strength from PSR J1852+0040. A con- 
vincing determination of the magnetic field strengths of 
the Cas A CCO and others will have to await discovery 
of their pulsations. 

4.4. Application to SN 1987 A 

It has long been recognized that the non-detection of 
the expected pulsar in the remnant of SN 1987A can be 
explained if the NS was born with a weak magnetic field 
or a long rotation period (Ogelman & Alpar 2004), and 
that several CCOs, including the one in Cas A, would be 
undetected if placed in SN 1987A (Graves et al. 2005). 
However, it was not established until now that a CCO 
can have essentially the same spin-down luminosity at 
birth that it has at an age of 400 or 10 4 yr, i.e., that it is 
born with both a weak B-field and a long P. This new re- 
sult means that a common type of neutron star can emit 
less than the observed limits from SN 1987A even if 100% 
of its spin-down power is reprocessed into IR emission by 
dust in the surrounding SN ejecta. The relevant luminos- 
ity limits to be satisfied for a point source inside the ring 
ofSN 1987A are L x (2-10keV) < 1.5 xlO 34 ergs s^ 1 cor- 
rected for extinction (Park et al. 2004), £(2900-9650 A) 
< 8 x 10 33 ergs s _1 corrected for dust absorption (Graves 
et al. 2005), and £(10 /xm) < 1.5 x 10 36 ergs s" 1 for dust 
emitting at T « 90 - 100 K (Bouchct et al. 2004). This 
mid-IR luminosity can be accounted for by radioactive 
decay of Ti, and therefore represents a very conser- 
vative upper limit on the spin-down power of an em- 
bedded pulsar. At an age of 10 — 20 yr, a cooling NS 
need emit only w 3 x 10 34 ergs s _1 of soft X-rays at 
a temperature of 2.5 x 10 6 K (Yakovlev et al. 2002), 
some of which is absorbed by SN ejecta or ISM in the 



LMC. The upper limit on spin-down power from PSR 
J1852+0040, E < 7 x 10 33 ergs s" 1 , is a challenge to 
detect in SN 1987A, and E < 1.5 x 10 32 ergs s" 1 from 
IE 1207.4-5209 would be impossible. So we conclude 
that a CCO is a promising model for the unseen NS in 
SN 1987A. 

5. CONCLUSIONS AND FUTURE WORK 

We obtained three new X-ray timing points on PSR 
J1852+0040 that, in combination with two previous ob- 
servations, reveal no significant change in spin period on 
time scales ranging from 1 week to 2.4 yr. In the dipolc 
spin-down formalism, this implies E < 7 x 10 33 ergs s _1 
and B p < 1.5 x 10 11 G. Such a low £?-field is not seen 
in young radio pulsars, and while it does not demand a 
special explanation, it does favor accretion as a source of 
the small, hot thermal region that is responsible for the 
highly modulated X-ray pulsations, which may exceed 
the spin-down power of the isolated NS, and strains any 
plausible interior cooling or external heating model. 

We showed that the observered spin period, 105 ms, in 
combination with a weak field, B p ~ 10 10 G, would al- 
low disk accretion in the propeller regime, which may 
be a more satisfactory model overall. If B p is even 
weaker, < 10 9 G, PSR J1852+0040 could be a "slow" 
accrctor. Considering also the properties of the CCO 
IE 1207.4—5209, we speculated that this class may have 
weak magnetic fields as a result of slow natal rotation, 
which enables them to accrete. While we do not deci- 
sively favor accretion over a rotation-powered pulsar, we 
emphasize that for PSR J1852+0040 a weak magnetic 
field is an unavoidable feature of either interpretation. 
Such low-i? births may be very common if they also oc- 
curred in Cas A and/or SN 1987A. 

A radio pulsar detection of PSR J1852+0040 would 
demonstrate cleanly that it is rotation powered, but our 
negative results so far are inconclusive because there are 
several reasons why a rotation-powered pulsar may not 
be detected in radio. Continued X-ray monitoring will 
provide a more sensitive test for spin-down or accretion 
torques, perhaps eventually in a definitive manner, and 
more evidence of whether or not the X-ray flux and pulse 
profile are variable. It is now important to obtain a 
phase-connected timing solution, the most efficient and 
sensitive method of measuring a small P. If obtained 
soon, a well-timed series of observations should allow a 
coherent solution to be fitted retroactively back to 2004. 
Also, considering the possibility that a weak magnetic 
field and accretion may be responsible for the absorption 
lines in the soft X-ray spectrum of IE 1207.4—5209, a 
deeper X-ray spectral study of PSR J1852+0040 may be 
revealing. Finally, in order to test the general applicabil- 
ity of these ideas to the class of CCOs, more sensitives 
searches for their pulsations are required. 

This investigation is based on observations obtained 
with XMM-Newton, an ESA science mission with instru- 
ments and contributions directly funded by ESA Member 
States and NASA. Support for this work was provided by 
NASA through XMM grant NNX06AH95G and Chandra 
Award SAO GO6-7048X issued by the Chandra X-ray 
Observatory Center, which is operated by the Smith- 
sonian Astrophysical Observatory for and on behalf of 
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NASA under contract NAS8-03060. The GBT is op- under cooperative agreement by Associated Universities, 
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